Asymptotic dynamics of short-waves in nonlinear dispersive models 
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The multiple-scale perturbation theory, well known for long-waves, is extended to the study of 
the far-field behaviour of short-waves, commonly called ripples. It is proved that the Benjamin- 
Bona-Mahony-Peregrine equation can propagates short-waves. This result contradict the Benjamin 
hypothesis that short-waves tends not to propagate in this model and close a part of the old con- 
troversy between Korteweg-de Vries and Benjamin-Bona-Mahony-Peregrine equations. We shown 
that a nonlinear (quadratic) Klein-Gordon type equation substitutes in a short-wave analysis the 
ubiquitous Korteweg-de Vries equation of long-wave approach. Moreover the kink solutions of 4 
and sine-Gordon equations are understood as an all orders asymptotic behaviour of short-waves. It 
is proved that the antikink solution of (f> 4 model which was never obtained perturbatively can be 
obtained by perturbation expansion in the wave-number k in the short-wave limit. 



Introduction. The method of multiples scales (or re- 
ductive perturbation method) is a very powerfull method 
to study a large number of physical phenomena, in par- 
ticular wave motion in nonlinear dispersive systems. It 
is well known that the far field dynamics of a long-wave 
(LW) with a small amplitude in a nonlinear and disper- 
sive system can be reduced almost always to a reduced 
set of model equations such as Boussinesq, Korteweg-de 
Vries (KdV), Modified KdV, etc 

The purpose of this work is to look for the far field be- 
haviour of short-waves (SW) in some nonlinear and dis- 
persive systems. In order to see how SW propagate we 
reformulated the method of multiples scales. The sys- 
tems we will consider are three. Firstly a system com- 
ing from hydrodynamic: The Bcnjamin-Bona-Mahony- 
Peregrine (BBMP) equation. Secondly we will study SW 
in two important classical relativistic field theory models: 
The 4 model (</> 4 ) and the sine-Gordon (SG) equation. 

We will prove for BBMP (@) that SW can build up 
the same soliton solution as obtained from LW || , hence 
rising the question of the unicity of the soliton descrip- 
tion. We will prove that the antikink (or kink) solution 
of </> 4 model (|2|) which can not be obtained as a pertur- 
bative solution in the non-linearity paramenter A occurs 
as a perturbative solution in the wave-number k in the 
SW limit. Moreover the kink solution of sine- Gordon (|) 
does not enter the classical LW perturbation scheme. We 
will prove that it appears as a perturbative solution for 
SW. 

Short- Wave Approach. We are interested in the prob- 
lem of the asymptotic dynamics of SW in nonlinear and 
dispersive systems. To perform this study, all degrees 
of dispersion of the system are taken into account in a 
Taylor expansion of the linear dispersion relation uj(k) 
around a large value of the wave- number k. The asymp- 
totic dynamics of SW for t — > oo is considered via the 
introduction of an infinite number of slow time variables 
T i i r 3 1 r 5 1 • ■ ■ an d a fast space variable £ following the ex- 
tension theory of Sandri j(| . The use of this fast variable 
and an infinite series of slow time variables constitute the 
first key of the SW approach. 



The solution is expanded in the form of a power se- 
ries in a small parameter e proportional to the inverse of 
the wave-number k. The perturbative series solution is 
secular. It is regularized (uniform expansion), through 
a renormalization of the frequency. This results to the 
celebrated Stokes' hypothesis on frequency-amplitude de- 
pendence in water waves . 

The Stokes' hypothesis is actually the second key tool 
of our approach. Indeed, for LW asymptotic description, 
the KdV §§ or MKdV g§ hierarchies occult the need 
of this tool, as they naturally provide the correct series 
expansion of the frequency. 

Basic Models. Hence the problem is the asymptotic 
behaviour of a SW in Benjamin-Bona-Mahony-Peregrine 
(BBMP) |lT| and in two classical relativistic nonlinear 
models: (frond sine-Gordon (SG) JL2||: 



(BBMP) u t + u x - u xxt = 3(u 2 ) x , 

(0 4 ) <f>xx - 4>tt = m 2 (j) - Xtf) 3 , 
^ 3 



(SG) 



fit 



—=sin\(- 



(1) 
(2) 

(3) 



The above models have quite different intrinsic char- 
acteristics. First of all SG is an integrable models while 
BBMP and </> 4 are not. Second the linear dispersion re- 
lation cu(k) has a finite limit as k — > oo (SW limit) for 
BBMP , while it is unbounded for SG and 4 . Indeed we 
have 



W(BBMP) 



1 + fc 2 ' 

W (SG) = 



(to 2 



(4) 
(5) 



However the phase and group velocities are all bounded 
in the SW limit k — > oo, which is a central point in this 
approach as indeed this very property allows the three 
models to sustain short waves. Then we face the problem 
of the nonlinear propagation of a SW, which is the object 
of this work. In the following, both nonintegrable sys- 
tems (BBMP and will be displayed in details, while 
the integrable one (SG) will be only sketched. 
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The BBMP Model. Let us consider a SW in (1) char- 
acterized by k = fc e _1 with k ~ C(l) an d e< 1. The 
plane wave solution of the linear problem u = exp i{kx — 
uj(k)t} inspires a fast variable £ = e~ 1 x and infinitely 
many slow time variables T2 n +i = e 2n+1 t (n = 0, 1, 2, ...), 
by expanding u> in powers of e. 

We assume the expansion 



u = u Q + e 2 u 2 + e 4 w 4 + ... . 

and suppose the extension U2 n = U2 n (C, Ti, T3, • • •), n 
0, 1, • • • , [0,p|. Then, the operators 



d_ _ ld_ 

dx e d( ' 

9 9 3 d 5 d 
dt 8t\ 8t 3 <9t5 



(G) 

1 = 

(7) 
(8) 



allow us to study the behavior of a short-wave for large 
time. 

BBMP gives at orders e _1 , e, e 3 , the equations (writ- 
ten only up e 3 ) 

-Uo,(ti +uq- 3ul = 0, (9) 
Lu 2 = uo,ti + mo,CCt 3! ( 10 ) 

Lm 4 = U 2)T1 - "0,r 3 + + M 0,CCT5 + 3 ( u 2)o (11) 



where L is the linear operator associated with 
L(v) = -uccn + ~ 6(um )c- 



(12) 



The unique solution of (^J) in the form uo(i]) with 
?7 = fcoC-^iTi-W3r3-w 5 r 5 ..., going to zero for \(\ -> 00 
is 



1 ^ 1 
u = -^sech 77, = - — • 

2 4fco 



(13) 



The values ti^u^, ■ ■ ■, corrections to the principal fre- 
quency loi (Stokes' hypoyhesis) are still free but will be 
determined later by the non-secularity requirement. 
The equation ( p^| ) for u 2 then reads: 

Lu2 = {Auj^k 2 — lo\ — 12w3fcQsech 2 ?7}sech 2 f7 tanhry, (14) 

and its two first right hand side terms are resonant (sec- 
ular producing terms) because |L3[ 

L(sech 2 7ytanh?7) = 0. 

These secular terms are eliminated by choosing 



w 3 



4fc 2 



1 



4 2 fc 3 ' 



Hence, (|l4|) yields the solution U2(rj) — A^ 1 k^ 2 u (ri). 

The equation (^) for u±(ri) contains secular producing 
terms originated by the first four terms in the right hand 
side. They can be eliminated choosing oj§ — — 4 k^ . 



The solution is u^n) = 4~ 2 fc^~ 4 ito(?7). This procedure 
can be repeated at any higher order n = 0, 1, 2, .. and we 
obtain recursively 



u 2 n(v) 



1 



W2n+1 



4 „+l fc 2n+i' 



(15) 



Next, nicely enough, not only the perturbative series 
solution ([3]) can be summed to give 



Ak 2 



n=0 



4fc 2 



1 



u (t}). 



(16) 



but also, by using u>2n+i, its argument 77 in the laboratory 
coordinates, which results as 



n=0 y ' 



kt 



Ak 2 -l 



(17) 



Therefore, this SW perturbation technique finally leads 
to the solution 



u(x, t) = — 



2k 2 



1 -4k 2 



sech 2 [fc(a: 



t 



4fc 2 



(18) 



This very expression, solution of BBMP, was obtained in 
H as an asymptotic limit of a LW of small amplitude. 
Thus, for t — > 00, the nonlinear dynamics of a SW (with 
an order one amplitude) and that of a LW (with small 
amplitude) are indistinguishable in BBMP. The equation 

is a nonlinear Klein-Gordon equation which sustitutes in 
this SW approach the classical Korteweg-de Vries of the 
LW approach. 

The 4> A Model. The topological antikink type so- 
lution of 4> A will be obtained by perturbation expan- 
sion starting from the constant solutions $0 — im/y/X. 
Hence we seek a solution 4>(rj) such that <fi — > ^fm/V\ for 

T) — > ±00. 

For rj < 0, the function u = <j) — m/V~X, goes to zero 
for 77-^—00 and satifies 



= —2m 2 u — imVXu 2 — Au 3 



(19) 



For unidirectional propagation the convenient fast vari- 
able C and the slow variables T2 n +i are in this case: 
( = e~ x (x — t), T\ = et,r 3 = e 3 t, ■ ■ ■. Expanding u ac- 
cording to u = € 2 (uq + e 2 U2 + e 4 W4 + •••)> the resulting 
equations are (up to e 6 ), 

Lu = 0, (20) 

Lu 2 = -2u Xt 3 + Wo,2n - 3mVAtto, (21) 
Lu 4 = u 2: 2n - 2u 2: Cr 3 - 2 "o,Ctb + 2w ,tit 3 

-6mVXulu 2 - Au 3 . (22) 
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with L in this case being the linear Klein-Gordon opera- 
tor 



L(v) = 2u Cti + 2m 2 



(23) 



We choose for the solution u$ of ( p0| ) the form uq = 
Bexp2r) with r\ = fcoC — (4fco) _1 m 2 Ti + W3T3 + • • • and 
B a constant. All linear terms at the right hand side of 
equations for U2( n -i) are secular. They can be eliminated 
choosing appropiately u>2 n -i, namely 



W2n-1 



rn 



2n 



n\{\-n)\2 n k 2n - 1 ' 



Next the solutions read 



(24) 



(25) 



With these values of cJ2n— l the series for 77 can again 
be summed as 



r\ — kx 



t, 



and also the perturbative series for u only if we choose 
B = -^2fc~ 2 . It leads for cf> = u + mj\f\ to 



777. / 77% 

<t> = --j= tanh{fcx - y k 2 + — t- logfc}, (26) 

To get this expression it is necessary to use the Fourier 
representation (x < 0) 



^2(-l) n+1 S n exp(2nx) = tanhx, 
n=0 

where S n are the Neumann's numbers (Sq — 1, S n — 
2,Vn = 1,2,3,...). 

The above solution is the antikink solution of 4 
(with an initial shift logfc/fc), which has not been ob- 
tained previously within another perturbation scheme. 



The expression y k 2 + can be interpreted as a non- 
linear frequency which defines the nonlinear group 
velocity 



dk 



k' 2 



(27) 



It is remarkable that the Lorentz invariance of ( pq ) is 
precisely related to that particular velocity, indeed 



m , m a 

7x tailh{ Vl ( vT 



t 



= ) -logfc}. (28) 



Note that the case ?y > in the perturbative series 
would simply yield the solution 0(— 77). 



The SG Model. Finally in the case of the sine-Gordon 
model (||), for = e(0 o + e 2 2 + e 4 04 + ••■ ), with (f> 2n 
functions of rj = fcoC + (2ko) m ti + CJ3T3 + ... where 
C = e~ l (x — t), n = ei,T3 = e 3 t, we obtain (up to 
order e 4 ) 



L(<t>o) = 0, 
L(h) 
L(0 4 ) 



A 3 

-2^0,Cr 3 + 0O,2ti - gj^O) 
'202, Cr 3 + 02, 2n - 20o, f T3 

'3! W2 + ^5!' 



20o, 



with L being in this case the operator 



L(v) = 2v, 



m 2 v. 



(29) 
(30) 

(31) 
(32) 



We choose for the solution of (|3^) the expression 0o = 
Cexp?7 with C a constant. As in the previous case all 
the linear terms at the right hand side of equations for 
02(n-i) ar e secular. They can be eliminated choosing 
W 2 (n-1) as 



W2n-1 = (-1) 



n+1 



,2;> 



n!(i 



Hence the solutions 



>2(n-l) r 



2 

ead 



(33) 



°2(n-l) 



-A^ I1 _ 1 C 2 "- 1 exp(2n- 1)77 



16' m 2 ("- 1 )(2n - 1, 
The series for summs for C = 4m / V~\ko and yields 
x _ 4m y^_^„exp[(2n + !)(??- logfc)] 



(34) 



n=0 



4m 



2n+ 1 



arctan |exp(fcx — \/k 2 — m 2 t — log fc) j . (35) 



In this case the Lorentz invariant form of ( |35| ) ap- 
pears as a function of the nonlinear phase velocity v — 
yjl — m 2 /k 2 as 



4m 

= —= arctan ^ exp 
V A 



:(ai — ui) — logfc 



(36) 



Conclusion and Comments We have applied a 
multiple-time version of the reductive perturbation 
method to study the solitary-wave and the kink-wave 
solutions of some nonlinear dispersive models. These so- 
lutions has already been known before. The alternative 
way gives here to obtain them shows that they repre- 
sented short-waves asymptotic dynamics (t — > 00) 

1 - BBMP for long-waves serves about the same pur- 
pose as KdV, whereas their behaviors in propagating 
short-waves can be expected to be rather different from 
KdV: from a linear analysis, BBMP does not propagate 
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short- waves while KdV amplify them Then our re- 
sult answers this old controversy on the relative relevance 
between KdV and BBMP. Q . Actually we proved that 
short-waves do propagate nonlinearly in BBMP models, 
and build up soliton-like solutions as t — + oo. 

2 - The antikink (or kink) solution of </> 4 model which 
cannot be obtained as a perturbative solution in A [ [l2|| , 
appears as a perturbative solution in k in the short-wave 
limit. 

3 - The eq. ( ]35| ) shows that the kink solution of SG is 
obtainable only from a short-wave dynammics, as indeed, 
the limit k — > gives rise to an imaginary argument. 

4 - An initial profile generically contains short-wave 
components which are usually neglected in favor of the 
long-wave components, which occurs as well through nu- 
merical discretisation of the models. As we have shown 
that the short-wave components asymptotically build up 
soliton solutions, the common understanding of a soliton 
as originating from long-wave is to be questionned. 

5 - It is worth noting finally that, within the long- 
wave approach, the nonlinear character of the solution is 
present already at the first orders, as indeed one usually 
finds the Boussinesq, KdV, MKdV, etc equations. This 
is not the case with the short-wave approach where usu- 
ally all orders (i.e. all times) are necessary to unveil the 
nonlinear character of the solution. 
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